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Abstract 

We present a closed bouncing universe model where the value of coupling constants is set by the dynamics of a ghost-like 
dilatonic scalar field. We show that adding a periodic potential for the scalar field leads to a cyclic Friedmann universe 
where the values of the couplings vary randomly from one cycle to the next. While the shuffling of values for the couplings 
happens during the bounce, within each cycle their time-dependence remains safely within present observational bounds 
for physically-motivated values of the model parameters. Our model presents an alternative to solutions of the fine 
tuning problem based on string landscape scenarios. 
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Introduction. A fundamental problem in particle physics 
and cosmology concerns the specification of the constants 
of nature, in particular the 19 free parameters of the Stan¬ 
dard Model. It appears that these parameters are fine- 
tuned to allow for the formation of complex structure and 
eventually life [T] . While the coupling constants of our uni¬ 
verse are not the only ones which could lead to such struc¬ 
tures, only some subset of all possible coupling constants 
could do so. Possible solutions require new physics at high 
energies, as is the case with superstring theory [5] . For ex¬ 
ample, the Heterotic string gives rise to a four dimensional 
chiral gauge theory with many of the ingredients to real¬ 
ize the Standard Model. However, these four dimensional 
compactifications present a landscape of vacua and cou¬ 
pling constants. The dynamics of strings in the early uni¬ 
verse were investigated in order to build models of string 
cosmology mil]. While it was the hope that string the¬ 
ory would univocally determine the measured couplings 
of the Standard Model, another approach emerged: the 
multiverse hypothesis [ills]. 

Eternal inflation generically predicts that while infla¬ 
tion ended in our local Hubble radius, it continues in other 
regions, triggering the emergence of a plethora of causally- 
disconnected bubble universes. If each bubble universe is 
endowed with different coupling constants-as generically 
realized in string theory-then one can use anthropic rea¬ 
soning to justify the values found within our cosmic hori¬ 
zon, given that we are here to ask the question. This 
marriage between eternal inflation and the landscape of 
possible perturbative string compactifications provides a 
resolution to the pressing question of fine tuning in mod¬ 
ern physics. One can, however, wonder whether there are 
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alternatives to the string landscape as a dynamical mech¬ 
anism to determine the couplings of the Standard Model. 

In this work, we propose a model to explain the ap¬ 
parent fine-tuning of coupling constants without recourse 
to the multiverse. We show that in a cyclic universe the 
fundamental constants can change pseudo-randomly from 
cycle to cycle. (We will qualify “pseudo” later.) Our cur¬ 
rent universe is then just the cycle which happens to con¬ 
tain a set of constants conducive to life. Cyclic universe 
models have previously been investigated as alternatives 
to inflation ■ The idea that different string vacua could 
be explored in different cycles has been suggested in the 
context of explaining the value of the cosmological con¬ 
stant [5]. A recent development in the path towards well- 
behaved cyclic cosmologies is the proposal of the anamor- 
phic universe |S]. This approach solves the problem of 
anisotropic instabilities which often plague bouncing mod¬ 
els. It also provides a mechanism for producing a nearly 
scale-invariant spectrum of perturbations. 

Here we will present a toy model for how a cyclic uni¬ 
verse with pseudo-randomly changing constants might be 
realized. One key ingredient is to promote all coupling 
constants to moduli fields, and dynamically demonstrate 
two features: i. During each bounce the coupling constants 
vary pseudo-randomly; ii. During the expansion phase in 
each cycle the time variation of the coupling constants re¬ 
main consistent with current observational bounds. For 
simplicity, we will focus on the gauge sector of the Stan¬ 
dard Model and propose how to generalize to the Yukawa 
sector in the conclusion. 

The Model. The possibility of a cyclic universe with chang¬ 
ing constants has been investigated before [10]. In that 
work, the bounce is caused by a free ghost scalar field 
whose kinetic energy is negative and scales as where 
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a{t) is the FRW scale factor. The ghost dilaton field de¬ 
termines the value of a coupling constant, in this case the 
electromagnetic coupling constant. The universe is also 
assumed to be closed and to contain radiation. These in¬ 
gredients allow for a series of closed universes separated 
by bounces. The value of the ghost field (and thus of the 
coupling) increases quickly and by the same amount dur¬ 
ing each bounce and then remains approximately constant 
during the following expansion/contraction cycle. The 
monotonically increasing coupling limits the feasibility of 
the model as a solution to the fine tuning problem. We 
note that while ghost fields remain problematic, we adopt 
the same phenomenological semi-classical approach as the 
authors in m , which is to avoid its quantization. In¬ 
deed, ghost fields have found widespread applications in 
field theory and cosmology, for example as candidates for 
phantom dark energy [TT] and k-essence inflation [H] . Ad¬ 
ditionally, in the anamorphic universe approach mentioned 
in the introduction, a kinetic term with the wrong sign can 
be rendered ghost free in the presence of a non-minimal 
coupling to gravity [9]. We are currently investigating 
whether our model can be embedded in the anamorphic 
framework and plan to report on this in future work. 

Our model incorporates a potential for the ghost field 
in a Friedmann universe. The action is 


S = J d'^x 

with 


- Y + 2V (V')] + Sgf 




( 1 ) 

( 2 ) 


where the coupling field for the i-th sector of the Stan¬ 
dard Model is gyj^^ = , with g^ constant, and M* 

some mass scale, which from here on we will take to be 
the Planck scale Mp. For clarity, we will focus on only one 
gauge sector; our approach is easily generalized to other 
sectors. With our metric signature, -I-,-I-), e = -|-1 

corresponds to a regular scalar field, while e = — 1 cor¬ 
responds to a ghost field. We take the potential to be 
periodic but negative, 


=-A‘^(l-bcos(V'//)). (3) 


The negativity of the potential ensures that there is no 
net cosmological constant during an expansion cycle, given 
that the negative kinetic energy density will drive the field 
to the potential maximum, where V (■0) = 0. The energy 
density and pressure of the field tp are 

A^(l-bcos(V’//)) (4) 

+ cos(?/'//)) (5) 

where / sets the energy scale as in axion-like models. 

The equation of motion for ip in an FRW spacetime is 


ip + 3Hip 


7 


sin(V'//) = 0, 


where H = df a. We assume that other relativistic degrees 
of freedom are modeled by a generic radiation term, so 
that the Friedmann equations are 

^ (_ 1^2 _ ^4(1 ^ - I’ 

(7) 

“ =(-7 + A^(l+ cos(V'//)) + ^) , (8) 

where Pro is the radiation energy density aX a = 1, K = 
±1,0 gives the spatial curvature and we have taken e = — 1. 

The hope is that the field ip will climb onto one of 
the potential maxima as the universe expands so the cou¬ 
pling constant that it determines will not change signif¬ 
icantly. As the universe contracts, the ip field acceler¬ 
ates. Its negative kinetic energy increases until it coun¬ 
teracts the radiation energy density and causes a bounce. 
At the bounce, the field is traveling quickly and can run 
across many maxima of the potential in both directions, 
resembling a sphaleron solution in electroweak baryogen- 
esis. The precise location in the potential where it settles 
will set up new initial conditions for the next bounce. The 
field will then traverse a different number of extrema the 
next time there is a bounce, possibly leading to a random 
walk among maxima over many cycles. (Our model can 
evade the Tolman problem that plague cyclic universes by 
adding interaction terms that create entropy via the mech¬ 
anism discovered in M) 

We will work in conformal time as the bounces occur 
over a longer period of conformal time than cosmic time 
making numerical solution easier. Writing Eqs. (§ and 
(|^ in dimensionless form in terms of conformal time we 
have 


= -2H^' + ^ sin(^'//); 

II IT I 

« =-^ + -7(1±cos(4'//)), 

a 3a 6 d 


(9) 

( 10 ) 


where = ip/Mp, H = a'ja, j3 = A^/p^o, / = f/Mp, 
and the dimensionless conformal time is ?) = {y/pPo/Mp)r], 
with primes denoting derivatives by fj and Mp = I/VSttG. 
The first Friedmann equation becomes 


^ -"“e T 


(l±cos(T//)) 


± 


3a2 


KM^ 

( 11 ) 
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When /? = 0, these equations reduce to the model of Bar- 
row et al. [To] and we have exact solutions 


, Va 


= 


9 ’ 
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(p) = X 1 ± Vl - 6A sin(77 ± go) 
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( 12 ) 

(13) 
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for constants A and 770 depending on initial conditions. The 
normalization of a is fixed by choosing the dimensionless 














curvature, KM^/p^o = +1. The maximum and minimum 


values of a are 


a 


max,min 


l(i±yTV15). 


(14) 


When ^ = 0 we can expand the solution about the bounce 
as 


1 ( 77 ) =. 


1 + 


1 


11 


(15) 


2 \ T^bounce / I 
with the bounce occurring at 77 = 0. We can plug this into 


Eq. (10) and set 77 = 0 to get 


riho 


= 


1 - 60 ^ • 


(16) 


This is a useful quantity since the timescale of the bounce 
determines how short the time steps of a numerical solver 
need to be in order to correctly go through the bounce. 
We would therefore also like to know this quantity when 
P 0. Assuming that something like the solution in Eq. 


(121 holds even when we include the potential, 'k moves 


quickly through field space at the bounce since a is small. 
The sinusoidal term in Eq. 0 therefore averages to zero 
during the bounce and we get back the equation of motion 
with no potential, whose solution is indeed given by Eq. 
(121. The cosine term in Eq. © also averages to zero and 
by setting "H = 0 we get an equation for the scale factor 
at the bounce. 


2/3amin + Sainin “ + A — 0. 


(17) 


Note that since the solution in Eq. (12) is now only valid 


near the bounce, the constant A is not as easily determined 
from initial conditions as it is for the case with no poten¬ 
tial. We can, however, use this to determine the bounce 


time since the cosine term in equation (101 also averages 


to zero and we can plug in the ansatz of Eq. (15) to get 


Vho 


— 


1 - - 3/3a: 


4 ’ 
min 


(18) 


where we have eliminated A. We can use this to check that 
for a given bounce, we are using a time step small enough 
to correctly capture the behavior. 

Equations 0 and ( [To| ) contain only two independent 
parameters, (3 and /. In exploring the space of solutions 
we should also consider different initial conditions. Since 
the equations are nonlinear, the dependence of the solu¬ 
tions on the initial conditions will be nontrivial. There 
are, in principle, four initial conditions to set, a(0), a'(0), 
4>(0) and 4>'(0). However, a(0) can be fixed using the first 
Friedmann equation (Eq. ©), and using the other ini¬ 
tial conditions. For simplicity, we start solutions at the 
maximum scale factor so a'(0) = 0 and set 47(0) = tt/ so 
that the potential energy vanishes initially. Then the only 
initial condition left to vary is 4''(0). 



t 


Figure 1: The behavior of the field plotted against cosmic time, 
t (in units Mp/y/pro). Parameters for this solution are /3 = 1, f = 
10“^, ^^(0) = 0.3, with a time step in dimensionless conformal time 
r; of 5 X 10~^. Inset: Enlargement of the behavior of ^ between two 
bounces. 


We solve equations 0 and ( [Io| ) numerically. We focus 
primarily on the behaviour of the field 4^ since this deter¬ 
mines the coupling constant in our model. We are looking 
for solutions with three main properties: 


1. 47 remains approximately constant during the expan¬ 
sion and contraction phases. 

2. 47 changes relatively quickly during the bounce phase. 

3. The change in 4t can change sign from bounce to 
bounce in a pseudo-random way. 

These three properties allow physical constants to be ap¬ 
proximately fixed during each cycle, but to undergo a 
pseudo-random walk over many cycles. 

We find that in order to obtain solutions with the de¬ 
sired properties, we should have /3 ~ 1 and / 1, or in 

terms of dimensionful quantities, ~ / <?; Mp. 

For these small values of /, the field 4^ crosses through 
many potential maxima during a bounce. Since the direc¬ 
tion of the change in 4* during the next bounce depends 
sensitively on where in the potential the field ends up af¬ 
ter the current bounce, the exact evolution becomes very 
sensitive to the time step used: small errors can build to 
the point where they change the direction of a jump in 
4> which changes the subsequent evolution substantially. 
However, this sensitivity to the time step only affects the 
precise sequence of jumps and not the general behavior. 
As long as the time step is chosen small enough compared 
to the bounce times (Eq. ([T^), the numerical solution will 
at least be representative of the true solution. 

In Figurej^we show an illustrative solution. The sharp 
changes in the field 4^ correspond to bounces, while the pe¬ 
riods where the field is comparatively constant correspond 
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to the expansion and contraction phases. The field ap¬ 
pears to undergo a pseudo-random walk and since it sets 
the coupling constant of the gauge field, gvMi, the pa¬ 
rameter space of coupling constants is explored over many 
cycles. In the inset we show the scale of oscillations in d/ 
during the expansion and contraction phases. We see that, 
as we would expect, the oscillations are of order / = 10“^, 
while the changes during the bounce have a magnitude 
of about 7. In fact, the magnitude of the changes during 
the bounce are basically independent of any parameters as 
long as the scale factor at the bounce is small. The scale 
of the oscillations away from the bounce, however, is given 
by the parameter /. The variation of coupling constants 
during the expansion and contraction phases can therefore 
be made arbitrarily small by choosing / sufficiently small. 

We can relate the change in tp during the expansion and 
contraction to the change in the coupling constant gyM 
during this time. The coupling varies as gyM = goe^^^^ = 
goe^. If varies on the order of and ^ I then 
the fractional change in the coupling constant will be of 
the order 


^9ym ^ ^ ^ 

gvM Mp 


(19) 


Since the parameter / sets the variation of tp away from 
the bounce, the fractional variation of gyM will be of the 
order f /Mp. Observations by Webb et al. suggest that the 
fine structure constant may have varied by 


-=-0.72 ±0.18 X 10"® (20) 

a 
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Figure 2: Autocorrelation as a function of lag k for three values of / 
as defined in equation \22\ . From top to bottom the values of / are 
0.01, 0.1 and 0.3. The circles are calculated from the model, while 
the crosses are for a true random walk for comparison. The y axis 
limits are the same for the top and bottom plots as for the middle 
plot. 


We then take the set of differences 

A, = «',+! - (21) 

and define the autocorrelation of the differences as 


since the early universe US]. In our model this would re¬ 
quire f /Mp ^ 10"®. 

During a bounce, since the value of ip changes by ap¬ 
proximately 7Mp, the coupling constant changes by a fac¬ 
tor of « 10^®. If the bare coupling go is of order one, 
then the gauge field would often become strongly coupled 
and could even become very strongly coupled, complicat¬ 
ing its dynamics. However, the bare coupling may very 
well be many orders of magnitude smaller than one, so 
that even with a large change of value during a bounce, 
the effective theory remains safely perturbative. Given 
the general approach of our proposal, a viable universe-in 
the sense of being able to produce astrophysical structures 
conducive to the emergence of life-would be one where the 
couplings remain safely within the perturbative regime so 
as to emulate the Standard Model. Either way, one can 
assume that the majority of the contribution to the rela¬ 
tivistic degrees of freedom remains in a thermal state, such 
that the energy density in radiation evolves smoothly from 
cycle to cycle. 

We would like to characterize the extent to which our 
solutions for the field are well modelled by a random 
walk. One way to do this is to calculate the autocorrelation 
between differences in 'k from one cycle to the next. As 
a representative value of 'k from each cycle we take the 
value when the scale factor reaches a maximum; call this 


Rk 


2 -^ 1—1 


( 22 ) 


For a true random walk this will always be small for k ^ 0. 
We plot the autocorrelation as a function of the lag k for 
three values of / in Figure and compare it that for a 
random walk. Clearly when / = 0.3Mp the values are 
correlated (bottom plot). In fact, the solution for tk is pe¬ 
riodic. When / = O.lMp, the autocorrelation is positive 
for k up to around ten. For / = O.OIMp (and smaller), the 
autocorrelation is indistinguishable from that of a random 
walk. This does not mean that the behavior is truly ran¬ 
dom; the dynamics are fundamentally deterministic. It 
does mean though, that a random walk is a good model 
for our solutions, and that its statistical properties will be 
similar. This justifies our use of the term pseudo-random 
walk. 

It is apparent that the parameter / is critical in deter¬ 
mining the behavior of solutions. The requirement that 
variations in the fine structure constant are small over 
the lifetime of the universe tells us that f /Mp ~ 10"® or 
smaller. Note that the regime where / ^ 1 is also where 
the varying coupling constant is well-modeled by a random 
walk. While we have considered only a single gauge sector, 
this approach can be generalized to multiple gauge fields 
with independent ghost fields ipi. If these fields are not 


4 








coupled (an interesting possibility), and each has a poten¬ 
tial width parameter, fi 1, then their respective cou¬ 
pling constants will undergo independent pseudo-random 
walks. As the universe progresses through many cycles, 
the coupling constants will explore the parameter space. 

Conclusion. In this work we have provided an alternative 
cosmological model to anthropic arguments in the string 
landscape scenario for explaining the values of the coupling 
constants of the Standard Model. Our toy model uses dila- 
ton fields which couple to the gauge sector of the Standard 
Model. We numerically demonstrated that while during 
the bounce the values of the coupling constants undergo 
a pseudo-random variation, they are stabilized during the 
expansion epoch of the universe. We showed that con¬ 
sistency with observations naturally favors randomness. 
Although the mechanism stands alone as an illustration 
of how to implement random changes in couplings in a 
bounce universe, it’s also motivated by string-theoretic re¬ 
alizations of the Standard Model where dilatons play the 
role of coupling constants in gauge sectors [TB]. A simi¬ 
lar procedure can be implemented for Yukawa couplings, 
promoting them to dilaton fields with periodic potentials. 
As with the gauge sector, we expect them to vary pseudo- 
randomly during the bounce, while remaining consistent 
with time-dependent observational bounds during the ex¬ 
pansion/contraction phases. We could loosely refer to this 
approach as a multiverse realized in time, as one consid¬ 
ers the variations of coupling constants over many expan¬ 
sion cycles. Within this framework, our cycle would be 
one where the couplings remain within the perturbative 
regime, emulating the Standard Model. 
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